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We review tlie literature on what classical physics has to say about the Meissner effect and 
the London equations. We discuss the relevance of the Bohr- van Leeuwen theorem for the 
perfect diamagnetism of superconductors and conclude that the theorem is based on invalid 
^ \ assumptions. We also point out results in the literature that show how magnetic flux ex- 

^ O . pulsion from a sample cooled to superconductivity can be understood as an approach to the 

' magnetostatic energy minimum. These results have been published several times but many text- 

books on magnetism still claim that there is no classical diamagnetism, and virtually all books on 
superconductivity repeat Meissner's 1933 statement that flux expulsion has no classical explanation. 
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• I. INTRODUCTION 

o , 

' It is now a century since superconductivity was discovered by Kammerlingh-Onnes in Leiden in 1911. From the 
^»-^' beginning, there was considerable interest from theoretical physicists. Unfortunately progress has been slow and one 

' can safely say that the phenomenon is still not completely understood, at least not in a fundamental reductionist 
sense. It is therefore important that the things that can be understood are correctly presented in textbooks. To 
the contrary, textbooks often repeat two myths which have become so ingrained in the minds of physicists that they 

. hamper progress. These are: 

QQ ' • There is no classical diamagnetism (Bohr 1911,^ van Leeuwen 1919^). 

[ • There is no classical explanation of flux expulsion from a superconductor (Meissner 1933^). 

I Although both statements have been disproved multiple times in the scientific literature, these myths continue to be 
spread. We hope this review will improve the situation. 

Meissner found it natural that a weak magnetic fleld could not penetrate a Type-I superconductor. The fact that 
this is already in conflict with classical physics according to the Bohr-van Leeuwen theorem is rarely mentioned. On 
the other hand, Meissner found it remarkable that a normal metal with a magnetic field inside will expel this fleld when 
cooled to superconductivity. This was claimed to have no classical explanation, and consequently superconductors 
k> , were not considered perfect conductors. 

5_j ■ We will not present any new results in this paper. Instead, we will review some of the strong evidence in the archival 
Ci ' literature demonstrating that the above two statements are false. We begin with the Bohr-van Leeuwen theorem and 

then turn to the classical explanation of flux expulsion. After discussing the arguments for the traditional point of 

view, we will demonstrate why they are wrong. 



II. THE BOHR- VAN LEEUWEN THEOREM AND CLASSICAL DIAMAGNETISM 



The Hamiltonian H, for a system of charged particles interacting via a (scalar) potential energy U is 

N 2 

F(r,,p,)=^^ + [/(r,.). (1) 
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The particles have masses ruj , position vectors rj , and momenta pj . The presence of an external magnetic field with 
vector potential A{r) alters the Hamiltonian to 

i=i 

where ej are the charges of masses mj . Using this equation one can show that statistical mechanics predicts that the 
energy — the thermal average of the Hamiltonian — does not depend on the external field. Hence the system exhibits 
neither a paramagnetic nor a diamagnetic response. 

In his 1911 doctoral dissertation, Niels Bohr^ used the above equation and concluded there is no magnetic response 
of a metal according to classical physics. In 1919 Hendrika Johanna van Leeuwen independently came to the same 
conclusion in her Leiden thesis. When her work was published in 1921, Miss van Leeuwen^ noted that similar 
conclusions had been reached by Bohr. 

Many books on magnetism refer to the above results as the Bohr-van Leeuwen theorem, or simply the van Leeuwen 
theorem. The theorem is often summarized as stating there is no classical magnetism. Since this is obvious in the case 
of spin or atomic angular momenta — the quantum phenomena responsible for paramagnetism — the more interesting 
conclusion is that classical statistical mechanics and electromagnetism cannot explain diamagnetism. A very thorough 
treatment of the Bohr-van Leeuwen theorem can be found in Van Vleck.— Other books, such as Mohn,— Getzlaff,— 
Aharonifi and Levy^^ mention the theorem more or less briefly. Weaknesses in classical derivations of diamagnetism 
in modern textbooks have been pointed out by O'Dell and Zia.^ A discussion of the theorem can also be found in 
the Feynman lectures-'^ (vol. II section 34-6), which points out that a constant external field will not do any work on 
a system of charges. Therefore, the energy of this system cannot depend on the external field. In addition, the fact 
that the diamagnetic response of most materials is very small lends empirical support to the theorem. 



A. The inadequacy of the basic assumptions 

If we accept the popular version of the Bohr- van Leeuwen theorem as true, then one must conclude that the perfect 
diamagnetism of superconductors cannot have a classical explanation. A study of the proof of the theorem, however, 
reveals that it is only valid under assumptions that do not necessarily hold. The assumed Hamiltonian of Eq. ^ 
only includes the vector potential of the external magnetic field. But it has been known since 1920 that the best 
Hamiltonian for a system of classical charged particles is the Darwin Hamiltonian.^^ In this Hamiltonian one takes 
into account the internal magnetic fields generated by the moving charged particles of the system itself, plus any 
corresponding interactions. 

The simplest way to see that the total energy of a system of charged particles depends on the external field is to 
note that the total energy includes a magnetic energy^ 

Eb = ^ I B\r)AV f{B, + B,fAV, (3) 



where Bg is the external magnetic field and is the internal field. Using the Biot-Savart law this field is given, to 
first order in v/c, by 

N N , , 

SjVj X {r-_rj) 
^■1 



B,ir) = E B,{r) = ^ f l.-..|3^ ■ (4) 



Equation ([3]) makes it obvious that to minimize the energy the internal field will be, as much as possible, equal in 
magnitude and opposite in direction to the external field. This gives diamagnetism. 

Charles Gallon Darwin was first to derive an approximate Lagrangian for a system of charged particles (neglecting 
radiation) that is correct to order (t^/c)^.— li^"— In a time independent external magnetic field, there is then a 



conserved Darwin energy given by 

N r 
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Uirk) + E,. (5) 



In this equation, Vj are velocity vectors, A^ is the external vector potential, E^ is the (constant) energy of the external 
magnetic field, and A^ is the internal vector potential given by 

AKr,-;r•.,..)^f:il -^■ + (^;^^^^^)^^^ (6) 
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where rkj = \rj — r^l and ekj = {rj — rk)/rkj- Although there is a corresponding Hamiltonian, it cannot be written 
in closed form. When the Darwin magnetic interactions are taken into account, the Bohr- van Leeuwen theorem is 
no longer valid because the magnetic fields of the moving charges will contribute to the total magnetic energy. The 
fact that this invalidates the Bohr-van Leeuwen theorem for superconductors was stated explicitly by Pfieiderer^^ in 
a letter to Nature in 1966. A more recent discussion of classical diamagnetism and the Darwin Hamiltonian is given 
by Essen.— 

B. What about Larmor's theorem? 

Larmor's theorem states that a spherically symmetric system of charged particles will start to rotate if an external 
magnetic field is turned on (see Landau and Lifshitai^ §45). The rotation of such a system produces a circulating 
current and thus a magnetic field. Simple calculations show that this field is of opposite direction to the external 
field and so the system is diamagnetic. This idea was first used by Langevin to derive diamagnetism. But from the 
point of view of the Bohr-van Leeuwen theorem this seems strange. In Feynman's lectures^^ the Bohr-van Leeuwen 
theorem is therefore not considered to be valid for systems that can rotate. As noted above, the problem is simply 
that the Bohr-van Leeuwen theorem does not take into account the magnetic field produced by the particles of the 
system itself. When this internal field is accounted for, diamagnetism follows naturally, whether for systems of atoms 
from Larmor's theoremi^ or for superconductors from the Darwin formalismi^ 

C. The Shanghai experiment - measuring a diamagnetic current? 

If it is true that the classical Hamiltonian for a system of charged particles predicts diamagnetism, then why is the 
phenomenon so weak and insignificant in most cases? Is there any evidence for classical diamagnetism for systems 
other than superconductors? As stressed by Mahajan,— plasmas are typically diamagnetic. However, plasmas are 
not usually in thermal equilibrium so it is difficult to reach any definite conclusions from them. An experiment on an 
electron gas in thermal equilibrium, performed by Xinyong Fu and Zitao Fu3^ in Shanghai, is therefore of considerable 
interest. 

Two electrodes of Ag-O-Cs side-by-side in a vacuum tube emit electrons at room temperature because of their 
low work function. If a magnetic field is imposed on this system, an asymmetry arises and electrons fiow from one 
electrode to the other. For a field strength of about 4 gauss a steady current of ~ 10""'^* A is measured at room 
temperature. The current grows with increasing field strength and changes direction as the polarity is reversed. 
The authors interpret this result as if the magnetic field acts as a Maxwell demon that can violate the second law of 
thermodynamics.^^ In view of statistical mechanics based on the Darwin Hamiltonian,^^ it is more natural to interpret 
this result as a diamagnetic response of the system. The current — ^just as the super-current of the Meissner effect — is 
due to a diamagnetic thermal equilibrium. This means that no useful work can be extracted from the system. 

III. ON THE ALLEGED INCONSISTENCY OF MAGNETIC FLUX EXPULSION WITH CLASSICAL 

PHYSICS 

Meissner and Ochsenfeld'^'^'^ discovered the Meissner effect in 1933. To their surprise a magnetic field was not only 
unable to penetrate a superconductor, it was also expelled from the interior of a conductor as it was cooled below 
its critical temperature. The first effect — ideal diamagnetism — seemed natural to them even though it violates the 
Bohr-van Leeuwen theorem. As already discussed, by violating this theorem ideal diamagnetism should have been 
considered a non-classical effect. The flux expulsion, on the other hand, was explicitly proclaimed by Meissner to 
have no classical explanation. Meissner does not give any arguments or references to support this statement, but 
according to Dahl^ the theoretical basis was Lippmann's theorem on the conservation of magnetic flux through an 
ideally conducting current loop. Later, Forrest^'^ and others have argued that the magnctohydrodynamic theorem on 
frozen-in flux lines also supports this notion. In this section we discuss these arguments and point out that they do 
not rule out a classical explanation of the Meissner effect. 

A. Lippmann's theorem 

Gabriel Lippmann (1845-1921), winner of the 1908 Nobel Prize in Physics, published a theorem in 1889 stating 
that the magnetic fiux through an ideally conducting current loop is conserved. In 1919, when superconductivity had 
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been discovered, Lippmann^S. again published this result in three different French journals (see Sauer-2S). Although 
the idea of flux conservation is considered highly fundamental^^'' nowadays references to Lippman's theorem are hard 
to find. But at the time Lippman's theorem was quite influential, and P. F. Dahl24 explains how this was one of the 
results that made the Meissner effect seem surprising — and non-classical — at the time of its discovery. 

The proof of Lippmann's theorem follows by noting that the self inductance L of a closed circuit, or loop, is related 
to the magnetic flux $ from the current in the loop via 

$ = cLq, (7) 

where c is the speed of light and q is the current through the circuit, an over dot denoting a time derivative (see 
Landau and Lifshita^^ vol. 8 §33). The equation of motion for a single loop electric circuit is 

Lq + Rq + C-^q = £{t), (8) 

where R is the resistance, C the capacitance, and £{t) the emf driving the current. If there is no resistance, no 
capacitance, and no emf, this equation becomes 

Lq = 0. (9) 
Therefore, if the self inductance L is constant, Eqs. (O and (O tell us that $ = constant. 



B. Lippmann's theorem and superconductors 

Although Lippmann's theorem is correct, its relevance for the prevention of flux expulsion is not clear. For super- 
conductors there are two points to consider, the assumption of zero emf, and the fact that constant flux does not 
imply constant magnetic field. We consider these points one at a time. 

Consider a superconducting sphere of radius r in a constant external magnetic field B^. When the sphere expels 
this field by generating (surface) currents that produce = ~Bc in its interior, the total magnetic energy is reduced. 
The magnetic energy change is 

A£e . -3 (_) ^, (10) 

or three times the initial interior magnetic energy.^ This energy is thus available for producing the emf required to 
generate currents in the sphere's interior. The assumption of Lippmann's theorem — that the emf is zero — is therefore 
not fulfilled. 

When a steady current flows through a fixed metal wire both the fiux and the magnetic field distribution are 
constant. On the other hand, when current flows in a loop in a conducting medium, a constant flux through the loop 
does not imply a constant magnetic field because the loop can change in size, shape, or location. Normally current 
loops are subject to forces that increase their radius (see e.g. Landau and LifshitZi^S, vol. 8 §34 Problem 4, or Essen^S 
Section 4.1). In view of this, Lippmann's theorem does not automatically imply that the magnetic field must be 
constant, even if the emf is zero. 

Other authors have reached similar conclusions. Mci and Liang^l carefully considered the electromagnetics of 
superconductors in 1991 and write, "Thus Meissner's experiment should be viewed through its time history instead 
of as a strictly dc event. In that case classical electromagnetic theory will be consistent with the Meissner effect." 



C. Frozen-in field lines 



The simplest derivation of the frozen-in field result for a conducting medium begins with Ohm's law j = <tE. If 
CT — > oo we must have E = to prevent infinite current. Faraday's law then gives dB/dt = — c(V x E) = which 
tells us the magnetic field B is constant (Forrest, Alfven and Falthammar'^^). 

When dealing with a conducting fluid, the equations of magnetohydrodynamics and the limit of infinite conductivity 
in Ohm's law gives rise to^ 

— = S/x{vxB). (11) 

This result tells us the magnetic field convects with the fiuid but does not dissipate. In addition, Eq. (jlip can be used 
to derive the following two statements that are often referred to as the frozen-in field theorem)^ 
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1. the magnetic flux through any closed curve moving with the fluid is constant, and 

2. a magnetic field line moving with the fluid remains a magnetic field line for all time. 
The first of these is a just a restatement of Lippmann's theorem for a current loop in a fiuid. 

D. Magnetic field lines in superconductors 

As shown in the previous section, Ohm's law with infinite conductivity implies there can be no change in the 
magnetic field since this would give rise to an infinite current density. It is, however, not physically meaningful to 
take the limit cr — >■ cx) in Ohm's law. In a medium of zero resistivity one must instead use the equation of motion for 
the charge carriers. One can then derive^SiM 

dj e^n „ e . „ , 

J7 = E+—3XB, (12) 

d r m mc 

for the time rate of change of the current density. The time derivative here is a convective, or material, time derivative. 
Ohm's law simply does not apply, and therefore it cannot imply that the magnetic field does not change. 

The second statement of the frozen-in field theorem can be questioned on the same grounds. Just as with Lippmann's 
theorem, the conclusion that the magnetic field remains constant does not follow, even if the statement is assumed 
valid. The magnetic field is determined by the density of magnetic field lines and constancy of this density requires 
that the conducting fluid is incompressible. One can easily imagine that the fiuid of superconducting electrons is 
compressible, and that magnetic pressure^^ pushes the fluid (with its magnetic field lines) to the surface of the 
material. Corroborating this point of view, Alfven and Falthammar— (section 5.4.2) state that "in low density 
plasmas the concept of frozen-in lines of force is questionable." 



E. A classical derivation of the London equations 

In 1981 W. Farrell Edwards^^ published a manuscript with the above title in Physical Review Letters. This caused 
an uproar of indignation and the journal later published three different criticisms of Edwards' workj^"— Incidentally, 
all of Edwards' critics restated (or implied) the textbook myth that the Meissner fiux expulsion does not have a 
classical explanation. In addition, the journal Nature published a study by John Bryan Taylor-2£ pointing out the 
faults of Edwards' derivation. Of course, Edwards is not the only one to publish an erroneous derivation of the London 
equations. A much earlier example is the derivation by Moored from 1976. The fact that various derivations have 
been wrong, of course, does not prove anything — as long as a correct derivation exists. 

In 1966 Nature published a classical explanation of the London equations by Pfleidererl^ which did not cause any 
comment. Indeed it has not been cited a single time, probably because it is very brief and cryptic. However, a classical 
derivation of the London equations and flux expulsion can be found in a classic textbook by the French Nobel laureate 
Pierre Gilles de Gennes.— Because this derivation should be more recognized, we repeat it here. 



F. de Gennes' derivation of flux expulsion 

The total energy of the relevant electrons in the superconductor is assumed to have three contributions: the 
condensation energy associated with the phase transition Eg^ the energy of the magnetic field Eb, and the kinetic 
energy of the moving superconducting electrons E^ ■ The total energy relevant to the problem is thus taken to be 

E = Es+EB+Ek. (13) 

The condensation energy is then assumed to be constant while the remaining two can vary in response to external 
field variations. The super-current density is written as 

j{r) ~ n{r) ev{r), (14) 

where n is the number density of superconducting electrons and v is their velocity, which gives a kinetic energy of 

Ek^ J ^nir)mv^r)dV - J ^^J^/WdF. (15) 
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By means of the Maxwell equation V x B = Airj/c and Eq. ([3]) for Eb, the total energy (|T3)) becomes 

E = Es + ^ J [B^ + A2(V X Bf] dV, (16) 

where we have assumed that n is constant in the region where there is current, and the London penetration depth is 
given by 

/ mc^ 

V Ane^n ^ ' 

Minimizing the energy in Eq. (|16p with respect to 15 then gives the London equation 

B + A^ V X (V X =0, (18) 

in one of its equivalent forms. Notice that this derivation utilizes no quantum concepts and does not contain Planck's 
constant. It is thus completely classical. A similar derivation has been published more recently by Badfa-Majos.^'^ 

The conclusion of de Gennes is clearly stated in his 1965 book^ (emphasis from the original): "The superconductor 
finds an equilibrium state where the sum of the kinetic and magnetic energies is minimum, and this state, for macro- 
scopic samples, corresponds to the expulsion of magnetic flux. " In spite of this, most textbooks continue to state that 
"flux expulsion has no classical explanation" as originally stated by Meissner and Ochsenfeld"^ and repeated in the 
influential monographs by Fritz London^ and Nobel laureate Max von Laue.^"^ As one textbook example, Ashcroft 
and Mermin^ explain that "perfect conductivity implies a time-independent magnetic field in the interior." 



G. A purely classical derivation from magnetostatics 



One might object that the electronic charge e is a microscopic constant, and that the London penetration depth A in 
most cases is so small that it seems to belong to the domain of microphysics. So even if quantum concepts do not enter 
explicitly, the above derivation does have microscopic elements. It is thus of interest that from a purely macroscopic 
point of view we can identify the kinetic energy of the conduction electrons solely with magnetic energy4^~ The 
energy that should be minimized would then be Eq. ([3]) 

Eb = ^J B^dV. (19) 

Unfortunately, this will not provide any information about the currents that are the sources of B. However, if we can 
neglect the contribution from fields at a sufficiently distant surface — i.e. when radiation is negligible — it is possible 
to rewrite this expression as 

E's = ^Jj-AdV, (20) 

where B — ^ x A. The idea is to then apply the variational principle to EjA — — Eb^ i.e. the expression 

dV. (21) 

Here the energy is written in terms of the field A and its source j . There are many ways of combining (|19|) and (|20|) 
to get an energy expression, but it turns out that (PT|) is the one that gives the simplest results. 

Starting from ([2T|) and adding the constraint V ■ j = 0, the integration is split into interior, surface, and exterior 
regions. The result is a theorem of classical magnetostatics that states, for a system of perfect conductors the magnetic 
field is zero in the interiors and all current flows on their surfaces (Fiolhais et al.--) This theorem is analogous to a 
similar result in electrostatics for electric fields and charge densities in conductors, sometimes referred to as Thomson's 
theorem (see Stratton,'*^ Sec. 2.11). Readers should consult the paper by Fiolhais et aL— for details. 

Many other authors have also reached the conclusion that the equilibrium state of a superconductor is in fact 
the state of minimum magnetic energy of an ideal conductor. Some of these are, in chronological order, Cullwickj— 
Pfleiderer,^^ Karlsson;^ Badia-Majos;^ Kudinov)^ and Mahajani^ For example, Cullwick writes "The well-known 
Meissner effect in pure superconductors is shown to be an expected rather than an unexpected phenomenon..." while 
Kudinov^ explains that it "is worth noting in this connection that the expulsion of a magnetic field to the periphery 
(the Meissner effect) also occurs in a classical collisionless gas of charged particles and that this happens solely because 
this state is energetically favorable." 



EjA = 



\j.A^l-iyxAf 
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IV. CONCLUSIONS 



The reader may get the impression from our investigations above that we consider superconductivity to be a 
classical phenomenon. Nothing could be further from the truth. As implied by the Ginzburg-Landau theory^ the 
BCS theory)^ and the Josephson effect^ the phenomenon is quantum mechanical to a large extent. Since quantum 
physics must lead to classical physics in some macroscopic limit, it must be possible to derive our classical result from 
a quantum perspective. Evans and Rickayzen^^ did indeed derive the equivalence of zero resistivity and the Meissner 
effect quantum mechanically, but did not discuss the classical limit. What we want to correct is the mis-statement 
that the Meissner effect proves that superconductors are "not just perfect conductors." According to basic physics 
and a large number of independent investigators, the specific phenomenon of flux expulsion follows naturally from 
classical physics and the zero resistance property of the superconductor — they are just perfect conductors. 

In conclusion, we have carefully examined the evidence for the oft repeated statements in textbooks that 1) there is 
no classical diamagnetism and 2) there is no classical explanation of magnetic flux expulsion. We have found that the 
theoretical arguments for these statements are not rigorous and that, for superconductors, these particular phenomena 
are in good agreement with the classical physics of ideal conductors. 
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